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^ ABSTRACT 
< 

We present Spitzer 70 and 160 /xm observations of the COSMOS Spitzer survey 

(S-COSMOS). The data processing techniques are discussed for the pubhcly released 
Q products consisting of images and source catalogs. We present accurate 70 and 160 /xm 

source counts of the COSMOS field and find reasonable agreement with measurements 
^ in other fields and with model predictions. The previously reported counts for GOODS- 

0-{ North and the extragalactic First Look Survey are updated with the latest calibration, 

Q and counts are measured based on the large area SWIRE survey to constrain the bright 

source counts. We measure an extragalactic confusion noise level of <Tc = 9-4 =b 3.3 mJy 
Ci (9 = 5) for the MIPS 160 /im band based on the deep S-COSMOS data and report an 

updated confusion noise level of cic = 0.35 ±0.15 mJy {q = 5) for the MIPS 70 /um band. 

> 

Subject headings: galaxies: evolution — galaxies: starburst — infrared: galaxies 
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1. Introduction 
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The Cosmic Evolution Survey (COSMOS) is a deep multi-wavelength wide-area (2 deg^) pro- 



gram for studying the evolution of galaxies and active galactic nuclei (AGN) (Scoville et al. 2007 1 . 
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The COSMOS Spitzer (S-COSMOS) survey is comprised of the Infrared Array Camera (IRAC, 3.6- 



8//m) and Multiband Imaging Photometer for Spitzer (MIPS, 24, 70, and 160 //m) data (Sanders et 
aL|2007[ ). This paper presents the far-infrared (FIR) 70 and 160 //m MIPS observations of the field. 
Although the mid-infrared (MIR) 24 fim array is more sensitive to the detection of distant galaxies 
than the MIPS-Germanium (MIPS-Ge) 70 and 160 nm detectors, the 24 ^m data are biased toward 
warm AGN and are affected by broad mid-infrared PAH emission and silicate absorption features 
redshifted into the band. The strong MIR spectral features along with the large variations of the 



FIR/MIR continuum ratios (e.g.. Dale et al. 20051 yield highly uncertain bolometric corrections. 



The long-wavelength 70 and 160 fim observations directly measure the FIR peak of the spectral- 
energy distributions (SEDs) for redshifts z < 1.5 and are key for constraining the total infrared 
luminosities and star-formation rates of galaxies within the COSMOS field. The MIPS 70 and 160 
^m data provide an important piece of the puzzle in the quest of understanding galaxy evolution. 

The goal of the data paper is to document the data products to facilitate the ongoing research 
of the COSMOS field. The data products and a description of the observations and data reduction 
are provided as part of the large public repository of multi-wavelength data for the COSMOS 
field {Hubble Space Telescope Advanced Camera for Surveys (ACS), Koekemoer et al. 2007; radio, 
Schinnerer et al. 2007; X-ray, Hasinger et al. 2007; and optical and near-infrared, Capak et al. 
2007). We present two scientific results here: (1) the 70 and 160 fim source counts, and (2) the 
measurement of the confusion noise for the MIPS 160 //m band. 



2. Observations 

The MIPS S-COSMOS observations were carried out in 4 campaigns from 2006 January 
through 2008 January (Table 1). The project {Spitzer programs 20070 and 30143) represents 
over 450 hr of MIPS observations. All observations were taken using the scan mapping mode dur- 
ing nominal "Cold" MIPS campaigns (telescope temperatures low enough to yield good quality 



160 nm. data). The initial observations taken in 2006 Cycle-2 are described in Sanders et al. (2007 1 . 
The Cycle-2 observations were comprised of a shallow wide-area (1.75 deg x 1.97 deg) survey to 
quantify the level of cirrus within the field and a small (0.5 deg x 0.33 deg) deep "test" field. After 
the successful completion of the Cycle-2 program, deep observations over the entire COSMOS field 
were carried out in Cycle-3 (depth of about 2800 s, 1350 s, and 270 s in the MIPS 24, 70, and 160 
bands respectively). 

The Cycle-3 astronomical observational requests (AORs) were optimized specifically for the 
MIPS-Ge bands, without compromising the 24 fim data. In contrast, some early MIPS programs of 
other groups were designed for the 24 /im band at the expense of the MIPS-Ge bands. The forward 
and return scan legs were offset by 148" which provides sufficient overlap for the 70 fim array. 
Cross-scan dither offsets of ±0, 42, 83" and in-scan dither offsets of ±0, 18" were used between 
multiple maps to account for the unusable parts of the MIPS-Ge arrays and the unobserved central 
row of the 160 /im footprint. The values of the cross-scan dithers were also chosen to avoid the 
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overlap of 24 /im readouts between consecutive maps. We carried out both forward and reverse 
scan maps to help characterize the long-term transients of the MIPS-Ge detectors. 

The majority of the data were taken in slow scan mode. At the slow scan speed (2.6"s~^), 
each AOR consisted of 4 scan legs of 1.5 deg. Each AOR mapped 1.5 degx592", and 10 AORs were 
used to map the field once. In total, 13 slow maps (29.5 hr per map) were carried out in Cycle-3, 
along with one map at the medium scan rate (12.5 hr) to complete the awarded time. Five AORs 
were lost due to satellite downlink issues in the second epoch of Cycle-3 and were re-observed in 
early 2008 (Cycle-3c, Table 1). 

The scan direction of MIPS is determined by the date of observation, and the observations 
were carried out on the days that minimized the zodiacal light. Since the field is near the ecliptic 
plane, the zodiacal background contributes significantly to the total noise budget for the MIPS 
24 and 70 /im bands. The zodiacal light is not significant at 160 /um. The galactic cirrus level is 



low in the direction of the COSMOS field (Sanders et al. 20071 and is not the dominant source of 



confusion noise within the MIPS-Ge bands. Figures 1 and 2 show the final MIPS 70 and 160 /im 
images, combining all epochs (Table 1). Since the 70 //m and 160 /im arrays are on the opposite 
side of the MIPS field of view, the overscan regions yield slightly non-symmetric coverage. For both 
Cycle-3a and Cycle-3b, the entire ACS field was observed, and the MIPS over-scan regions provide 
coverage for the IRAC data outside of the ACS field. 



Data Reduction 



The raw MIPS-Germanium 70 and 160 /im (MIPS-Ge) data were downloaded from the Spitzer 
Science Center (SSC) archive and were reduced using the Germanium Reprocessing Tools (GeRT, 
version 20060415) and additional specialized scripts developed for processing the MIPS-Ge survey 
data. The GeRT uses an offline version of the SSC pipeline to produce the basic calibrated data 



products (BCDs). The basic MIPS-Ge processing steps are discussed by Gordon et al. (20051 



and within the GeRT documentation. The processing for S-COSMOS made use of lessons learned 



from the processing of the extragalactic First Look Survey (xFLS, Frayer et al. 2006a I and the 



deep GOODS-North observations ( Frayer et al.|[2006b |. Additional processing enhancements were 
derived here using the S-COSMOS data. 

The main processing steps were carried out in the following order: (1) calculation of the 
data ramp slope, (2) stimulator-flash interpolation, (3) improved stimulator-flash response solution, 
(4) calibration of the slope image to yield the BCD product, (5) enhanced flltering of the BCD 
product, (6) data co-addition, (7) identiflcation of bright sources, (8) masking bright sources and 
re-calculation of the flltering corrections, (9) flnal data co-addition, and (10) flnal source extraction. 
For comparison, the on-line SSC pipeline only performs steps (1), (2), (4), basic flltering, and step 
(6), and the GeRT includes software for steps (1), (2), (4), and (8). The 10 processing steps are 
summarized in the following sub-sections. 



- 4 - 



3.1. Basic Processing 



The optimal processing for step (1) depends on the background, the length of the data ramp 
(MIPS-Ge data are recorded with non-destructive reads sampled at 0.131 s), and the rate of cosmic 
rays at the time of the observations. The SSC on-line pipeline is tuned for the short data ramps (3 s 
and 4 s). We tuned the pipeline modules (cosmic ray detection and removal and slope estimation) 
of step (1) to minimize the noise level for the longer 10 s ramps of the S-COSMOS data. The tuning 
solutions for 70 fiui are similar to those derived for the GOODS-North photometry data which have 



the same ramp length (Frayer et al. 2006b I . MIPS-Ge uses stimulator-flash observations to track the 



response of each detector as a function of time. After the calculation of the initial stimulator-flash 
solution, which is basically a linear interpolation between stimulator- flash measurements (step 2), 
we removed outlier values and re-derived a smoothed stimulator-flash solution (step 3). For 70 fim, 
the stimulator-flash response function was smoothed by about 2 minutes (slightly longer than the 
stimulator-flash cycle) to provide the lowest noise. For 160 fim, the solution was smoothed by 
about 8 minutes to yield the best results. 

In step 4, the BCD data are calibrated as BCD{t) = FC[U{t) / SR{t) - DARK]/IC, where 
U{t) is the uncalibrated slope image, SR{t) is the stimulator-flash response solution derived in step 
(3), and DARK is the dark calibration file. The IC calibration file is the illumination correction 



which corrects for the flat-field response and the non- uniformity of the stimulator flash (Gordon 



et al. 20051. The flux conversion factor {FC) converts the instrument units into physical surface 
brightness units of MJysr"^. For self-consistency, we adopt the same DARK and IC files used 



for the official calibration of the MIPS instrument (Engelbracht et al. 2007 Gordon et al 



Stansberry et al. 2007 1 . The cahbration FC values of 702 MJy sr~^ per MIPS-70-unit (Gordon et al 



2007 



2007[ ) and 41.7 MJysr"! per MIPS-160-unit ( [Stansberry et al.||2007| ) are adopted. The calibration 
of MIPS is based on stellar SEDs {S y (X u ). We have placed the data on a constant u Sy scale by 
dividing the data by the color correction factors of 0.918 and 0.959 for the 70 and 160 ^m bands. 



respectively (Stansberry et al. 2007 1 . These color corrections are appropriate (accurate to better 
than 2%) for a wide range of galaxy and AGN SEDs {Sy oc a = 0-3) across the filter bandpass. 



3.2. Data Filtering 

Optimization of the processing steps (1-4) provided sensitivity improvements of about 20% 
in comparison to the default parameters, while the filtering (step 5) can provide more than a 
factor of two improvement in point-source sensitivity. Filtering is key in the removal of systematic 
instrumental effects which impede the ability to integrate down with deep observations. The two 
main artifacts impacting MIPS-Ge data are the latents due to the stimulator fiashes and variations 
of the slow response (> 2min) as a function of time. The slow response is removed at 70 and 
160 /Lim by subtracting a running median per detector as a function of time, i.e., a high-pass filter. 
The latent artifacts due to the stimulator fiashes are not fully removed by a simple high-pass filter. 
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since these variations occur on shorter time-scales. At 70 fim the stimulator-flash artifacts are 
correlated by column. Since the scan direction is nearly along the columns of the array, these 
artifacts contribute to the streaking within the maps if not corrected. We remove the column 
residuals by subtracting the median of the values along each column for every BCD at 70 //m. The 
combination of the column median filter and a high-pass median time filter per detector removes 
the instrumental artifacts at 70 fj,m. 

There is no equivalent column filter at 160 /um to remove the high-frequency (short-time 
scale) latent images introduced by the stimulator-flashes. Fortunately at 160 fj,m, these artifacts 
are repeatable and can be determined by stacking the data as a function of DCENUM (Data- 
Collection-Event Number) within the stimulator-flash cycle. Since the scan-mirror position also 
varies with DCENUM within the stimulator-flash cycle, we stacked the data per AOR for each 
scan-mirror position and took the median value of the stack to derive the correction as a function 
of detector and scan-mirror position. These corrections were subtracted from the BCDs to remove 
artifacts due to the stimulator- flashes. The combination of this stacking correction and a high-pass 
median time filter removes the instrumental artifacts at 160 /xm. 

With the S-COSMOS data, we tuned the filtering techniques (step 5) to minimize the noise 
in the output maps. For both 70 and 160 ^m, the short-time scale stimulator-flash artifacts were 
removed first, followed by a high-pass median time filter to remove the longer time-scale transients. 
We adopted a time filter width of 12 frames (2 minutes) to yield the best results. After the 
initial filtering, the brightest sources have negative "side-lobes" in the maps since bright sources 
bias the calculation of the median for neighboring pixels. To remove these filtering artifacts, the 
filtering was done in two passes. The data from the first filtering pass (step 5) were co-added 
(step 6), and sources were extracted (step 7) to find the location of the bright sources. The source 
positions within the original BCDs were masked and new filtering corrections were calculated in 
a second pass, ignoring the pixels containing sources. This two-pass filtering technique minimizes 
the artifacts while preserving the point-source calibration. After the second filtering pass (step 8), 
the data were co-added to produce the final maps (step 9, Sec. 3.3) and sources were extracted to 
produce the catalogs (step 10, Sec. 3.4). 



3.3. Imaging 



The SSC mosaicking software (MOPEX, Version 16.3.7, Makovoz &; Marleau 20051 was used 
to combine the data and make the final images. A fast plane-to-plane coordinate transformation 



method was used to project the data onto the sky (Makovoz 2004 1 with the default MOPEX 



interpolation scheme. We carried out the imaging steps following the techniques discussed for the 
xFLS MIPS-Ge data ( Frayer et al.|200"6a |. An important improvement available after the processing 
of the xFLS data is a more robust outlier rejection technique within MOPEX. The updated method 
rejects data around the median of a data stack for each sky pixel instead of rejecting data with 
respect to the average of the data stack. This enables more aggressive outlier rejection without 
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compromising the calibration of point sources. The best sensitivity was obtained using rejection 
thresholds of ±2.5(7. The new outlier rejection method improved the sensitivity in the maps by 
about 5%. 

The pipeline uncertainties of the BCDs and integration times were not used as weights in the 
co-addition of the data; all of the data (155,411 BCDs in total) flagged as good were given equal 
weight. Bad data flagged on a detector basis during the pipeline processing or identified as outliers 
by MOPEX were not included. The vast majority of the data were taken with 10 s integrations 
and no correction is needed for the 4 s BCDs (all BCDs are calibrated correctly in MJys~^). The 
pipeline uncertainties (calculated from the formal error propagation of the pipeline steps) do not 
fully represent the actual noise characteristics of the data, and underestimate the real noise slightly 
in low background regions. The main utility of the pipeline uncertainties for these data is to provide 
a lower limit to the input noise for the MOPEX outlier rejection algorithm. 

3.4. Source Detection and Extraction 

Sources were detected and extracted from the final images using the Astronomical Point-Source 
Extraction (APEX) tools within the MOPEX software package and using additional specialized 
scripts. For optimal source detection and extraction, accurate background subtraction and noise 
estimates are needed. The filtering process (steps 5&;8) yields a small systematic negative offset of 
the background level in the image. This offset was estimated by taking the median of the image 
within the central regions after masking sources detected at levels greater than Sa. After the 
removal of the global offset level (—0.05 MJysr"^ at 70 /xm and —0.07 MJysr"^ at 160 /itm), there 
are still local background fluctuations across the image depending on the local density of sources 
(both from detected sources and un-detected/confused sources associated with infrared galaxies 
seen at 24 /im). The local background level was derived by taking a median within a box around 
each pixel (after masking sources detected at levels greater than 3a). For source detection, we used 
a small box with a linear size of 5 Full- Width Half-Maximum (FWHM) widths to remove the local 
background, and for source extraction (fitting) we used a larger background box with a size of 9-10 
FWHM widths to conserve the calibration. 

Several types of noise images can be produced by the MOPEX and APEX software, but none 
are optimal for these data. The "std" noise image produced by MOPEX represents the empirical 
scatter from the repeated observations per sky pixel divided by the square-root of the number of 
good observations. For deep observations the "std" file underestimates the true noise since it does 
not account for the pixel-to-pixel correlated noise or the confusion noise. APEX computes the 
spatial pixel-to-pixel noise ("noise" file) by calculating the noise within a local box surrounding 
each pixel after rejecting positive outliers. The output "noise" image has variations that depend on 
the outlier parameter and the number of sources within the local box. To avoid local biases due to 
sources, we derived a "noise" image after the extraction of sources, adopting a box size with a linear 
scale of 9-10 FWHM widths (the same size used for the local background subtraction for source 
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fitting) . To preserve both the small and large scale spatial variations of the uncertainty across the 
science image, both the "std" (representing small-scale variations) and "noise" (representing large- 
scale variations) files were used. The quality of "std" and "noise" images was first improved by 
smoothing the images by 1 and 3 FWHM widths, respectively, and then combined in quadrature. 
Equal weights were given to the "std" and "noise" images at 70 fxm. Since the confusion noise 
is important at 160 fim, the weights at 160 ^m were based on the relative contributions of the 
instrument and confusion noise (Sec. 5.2). The std image was given a weight corresponding to the 
instrument noise, and the noise image was given a weight corresponding to the confusion noise to 
produce the combined uncertainty file at 160 fim. The median level of the combined uncertainty 
file was scaled to match the total average noise level derived from fitting the la width of a Gaussian 
to the data histogram of the image after source extraction. Source extraction and noise estimates 
were repeated until the results converged. 

After proper background subtraction and deriving an accurate uncertainty image, sources 
were detected using the APEX peak algorithm. Peaks with a signal-to-noise ratio (SNR) of greater 
than 3 were fitted using the Point-source Response Function (PRE, which includes the effects of the 
detector size and the adopted sub-sampling of the detectors) image. At the spatial resolution of the 
MIPS-Ge bands, the PRE is stable. At 70 /^m we adopt the PRE (EWHM = 18.6") made previously 



from the xELS data (Erayer et al. 2006a). At 160 /xm we made a new PRE using the COSMOS, 
xELS, and the data from all of the fields from the Spitzer Wide-area Infrared Extragalactic Legacy 
Survey (SWIRE, [Lonsdale et al.||2004 |. None of the individual fields has a large number of isolated 



160 fxm sources with high SNR to make a high-quality PRE. At 160 /Ltm, the brightest sources (> 2 



Jy) cannot be used due to the non-linearity of the detectors (Stansberry et al. 2007 1 . In total, we 
used 149 isolated sources from the Spitzer 160 fim surveys which have SNR > 30 and S'160 < 1 JjQ 
to produce an empirical PRE (EWHM = 39"). 

The detection and extraction of sources was done in multiple passes. Initial source lists were 
visually analyzed and cleaned to remove sources with low coverage and potential spurious sources 
in the Airy rings around bright sources. The detection table was updated to remove bad sources, 
and source fitting was redone. In the few cases where the detection routine failed in the proper de- 
blending of two or more sources, the detection table was modified, and source fitting measurements 
were redone using better initial source positions. 



3.5. Astrometry 

The absolute pointing reconstruction of the Spitzer telescope is typically better than 1". The 
pointing uncertainties are much less than the large EWHM width of the data (18.6" at 70 fim and 



^Throughout this paper the flux densities of the MIPS bands are deflned as 5160 = 5^(155. 9^m), S70 = 
5'„(71.4/xm), and S24 = S„{23.7fj.m). 
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39" at 160 /um). We verified the pointing solutions by comparing the positions of the 70 sources 
with their VLA radio counterparts (Schinnerer et al. 2007). For the approximately 400 sources 
detected at SNR> 5 at both 1.4 GHz and 70 /im, we find an average positional difference {Spitzer 
- VLA) of ARA = -0.08 ± 0.08" and ADec = 0.17 ± 0.08". These smah offsets are well within the 
scatter measured for individual sources of 1-2", so no positional corrections were made to the data. 



3.6. Calibration 



MIPS is calibrated assuming the point-spread functions (PSF, which does not include the 
effect of the detector size) calculated using the Spitzer TinyTim (S-TinyTim) models (Krist 2002). 



The modeled PSFs have been shown to match the observations for all bands (Engelbracht et al. 



2007 Gordon et al. 2007 Stansberry et al. 2007 1 . An important aspect of the MIPS PSFs is that 



significant flux arises outside of the first Airy ring on large spatial scales. We corrected the derived 
flux densities for emission outside of the adopted PRF images using the S-TinyTim models. For the 
adopted PRFs (87" x 87" at 70 /xm and 190" x 190" at 160 /im), we derive an aperture correction of 
1.15 at both wavelengths based on modeled PSFs assuming a constant v Sy SED. Empirically, it is 
difficult to measure the level of emission outside of the ffist Airy ring for MIPS-Ge data. However, 
the calibration factors assume the modeled PSFs out to very large spatial scales, and we must 
apply this correction for consistency. We verified the calibration consistency of our techniques 
using archived calibration observations taken over the lifetime of the mission. For comparison, 
observations of the calibration star HD180711 (S'70 = 447.4 mJy, Gordon et al. 2007) were used 



at 70 //m, while observations of the ULIRG IRAS 03538-6432 were used at 160 /xm (Stansberry 



et al. 2007 Klaas et al. 2001, (5160 = 1.04 Jy)). Using an aperture correction for a modeled PSF 
with a stellar SED and similar reduction techniques carried out for S-COSMOS, we measure a flux 
density of 450 it 18mJy for HD180711. This is consistent with expectations (flux density ratio of 
1.01 lb 0.04). At 160 /im we measure a flux density ratio of 0.97 it 0.05 compared to expectations 
for IRAS 03538-6432 using similar reduction techniques carried out for S-COSMOS. These results 
suggest that the calibration of the S-COSMOS data is consistent with the MIPS calibration papers; 
the absolute calibration uncertainty of MIPS is 5% at 70 fim ( [Gordon et al][2007l ) and 12% at 160 
lim. ([Stansberry et al.|[2007l). 



3.7. Corrections for Completeness and Eddington Bias 

The completeness levels were estimated by simulations and calculations based on the SNR 
threshold. For the simulations, sources were injected at random locations into the images and 
extracted using the same techniques adopted for the production of the catalogs. For high SNR 
> 5.0, the SNR threshold itself is the dominant effect in the determination of the completeness 
level. At lower SNR, other effects such as the details associated with source detection and fitting 
become significant. The completeness levels vary significantly across the image as a function of 
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coverage and flux density. Based on the SNR > 5.0 threshold, the average effective completeness 
level across the image can be calculated as a function of flux density as the fractional area within 
the image for which S* + / > 5. Of/, where S is the flux density, / is the science image, and U is 
the uncertainty image. Figures 3&:4 show completeness calculations as a function of flux density 
and coverage. The simulations match the expected curves very well at 70 for both the deep 
and typical regions within the image. For the deep region at 160 /xm, the simulated completeness 
values are relatively noisy and are lower than expected for 60-90 mJy, potentially due to the eflFects 
of confusion and/or the small number of independent beams within the simulated area. The 
completeness calculations at 160 ^m for the nominal coverage of greater than 20 are similar to the 
simulations and the calculations for the typical range of coverage values (25-33). For the derivation 
of the S-COSMOS source counts (Sec. 5.1), we adopt the completeness curves for coverages greater 
than the nominal values of 100 and 20 at 70 and 160 /xm, respectively (solid lines, Fig. 3&;4). 

In addition to the completeness corrections, the counts are affected by the Eddington bias (flux 
boosting) . At faint flux densities the observed flux densities are slightly higher on average than the 
true flux densities since sources on positive noise features are preferentially selected. This effect of 
"flux boosting" is not as important for SNR > 5 as for lower SNR, but would still yield a small 
systematic bias in the measured counts if not corrected. The same simulations used to help quantify 
the completeness corrections were used to check the importance of the Eddington bias. The ratio 
of the observed to input flux densities (S'obs/'S'true) of each bin were measured for both the deep 
region and the wide area with "full-blown" simulations that injected sources at random positions 
within the image (one at a time) and exacted the output flux densities using the same methods 
adopted to produce the source catalogs. Given the variation of coverage in the data, it is difficult 
to obtain sufficient statistics as a function of ffux density and coverage across the image using 
these full-blown simulations (as was the case for the completeness estimates). Instead, we carried 
out simpler calculations that are consistent with the full-blown simulations, but are signiflcantly 
more accurate. Input flux densities (-Strue) with the same power-law distribution as the observed 
source counts were randomly added to the noise distribution given by the uncertainty image, and 
the output observed flux densities were derived. The effective Eddington bias of Sohs/Stmc was 
calculated over the entire image for each flux density bin for the adopted SNR > 5 cut and used 
to correct the observed flux densities for the derivation of the counts (Sec. 5.1). This method fully 
accounts for the variation of coverage and noise across the image. 



4. Products 

The S-COSMOS 70 and 160 fim products are available on-line at the NASA Infrared Science 
Archive (IRS A) at the Infrared Processing and Analysis Center (IPAC). The products described 

here are version 3 of the data. Version 1 and version 2 were early quick-look products based on 
simplified reductions of sub-sets of the observations. Version 3 represents the combination of all 
data from the S-COSMOS MIPS program and is the first version to be calibrated properly and 
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processed with the best known data reduction techniques. 

4.1. Images 

The science images have been background subtracted with the removal of a global offset, but 
the local background fluctuations have not been removed (Sec. 3.4). The images are in surface- 
brightness units of MJy sr~^ assuming the latest calibration and have been color-corrected to match 
SEDs with a constant v Su scale (Sec. 3.1). The uncertainty images have also been color-corrected 
and are in units of MJysr"^. The uncertainty images (la) represent both the small-scale and 
large-scale spatial noise properties associated with the science image (Sec. 3.4). 

The coverage maps provide the effective number of observations (after data editing) per point 
on the sky for the science images. At 70 /um for the nominal coverage of greater than 100, the 
median effective exposure time is 1350 s while the ultra-deep region with a coverage of greater than 
250 has an exposure time of 2800 s. At 160 /xm for the nominal coverage of greater than 20, the 
median effective exposure time is 273 s, while the ultra-deep region with a coverage of greater than 
50 has an effective exposure time of 567 s. 

Table 2 lists the image properties and sensitivities. The sensitivities represent the total noise, 
including confusion. The surface-brightness noise for the adopted pixel scale was derived from 
fitting the la width of a Gaussian to the data histogram of the image after source extraction. To 
derive the effective point-source noise, including the effects of correlated noise between the pixels, 
we carried out aperture measurements at random locations within the residual image after the 
extraction of sources. We derive conversion factors between the point-source noise and surface- 
brightness noise of 13.3 it 1.3mJy (MJysr~^)~^ and 74 it 7mJy (MJy sr^"'^)^^ for the 70 fim and 
160 ;uni images respectively. The average point-source noise (la) is 1.7mJy and 13mJy at 70 and 
160 //m, respectively. These point-source values include the aperture correction of 1.15 for emission 
outside of the measured PRF (Sec. 3.6). 

4.2. Catalogs 

Point-source catalogs were made for a SNR > 5.0. In total, 1512 sources at 70 /xm (Table 3) 
and 499 sources at 160 (Table 4) are cataloged. The catalogs are single-band source lists and 
arc independent from each other and the MIPS 24 /xm data. Although the catalogs are not biased 
by data at other wavelengths, we did use the 24 ^m and radio data to confirm the reliability of 
the catalogs. Within the central area cataloged at 24 /im, only seven 70 //m sources do not have 
a 24 fim counterpart (524 > GO/iJy) within the approximate 70 /xm beam radius of 9" (which 
corresponds to a relatively large chance positional coincident of about 50% for a 24 fim source) . Of 
these seven, four are associated with a radio source, two do not have a radio counterpart, and one 
is outside the radio coverage. For the three cases without a current 24 /xm or radio counterpart, the 



- 11 - 



70 fim position is located between a blend of two or three 24 fiui sources. These blended sources 
could represent valid detections at 70 /xm. All 160 fim sources have possible 24 //m counterparts, 
which is not unexpected given the low spatial resolution of the 160 fj,m data (where on average 
there are about three 24 /im sources per 160 fan beam). Although we find no obvious spurious 
detections within the catalogs, users should be cautious and check the images when comparing 
catalogs at different wavelengths due to potential source blending. Tables 3&:4 show the format for 
an example portion of the S-COSMOS 70 fim and 160 ^m catalogs published in the on-line edition 
of the Journal. 

PRF-fitting flux densities and aperture measurements were made using the APEX software. 
The aperture and PRF measurements are in statistical agreement. Since PRF measurements are 



significantly more accurate for faint sources (e.g., Frayer et al. 2006a I, the PRF flux densities are 
used for the vast majority of sources (flag of "p"). Aperture measurements are adopted for sources 
not reasonably well fitted by the PRF (extended sources or unresolved blends) and are given a flag 
of "a" within the catalog. Resolved blends fitted well by multiple PRFs are given a flag of "p". 
For the one blend of two resolved sources which is not fitted well by two point sources at 70 /im, 
the total aperture flux is divided among the two components based on their relative peaks (flag 
of "ap"). For consistency with the adopted calibration of MIPS, the measurements include the 
aperture correction of 1.15 for emission outside of the measured PRF (Sec. 3.6). No corrections 



for flux non-linearity have been made (Stansberry et al. 2007). Only one source (SCOSMOS160 
J100027. 0-1-032226, 5160 ^ 11 Jy) is bright enough to be significantly impacted by flux non-linearity 
(the only source at 160 /xm with an aperture measurement, flag of "a"); its flux density should be 
treated with caution. 

The errors on the fitted flux densities derived by APEX are not used since they are underes- 
timated by about a factor of 3 for these data. We adopt errors based on the SNR measurements 
which represent the fitted peak flux density divided by the uncertainty image at the location of 
the source. The errors on the measured flux densities {S) represent the random errors given by 
the SNR ratio combined with the systematic calibration uncertainty (ecai)- The flux density error 
a{S) = (1/SNR e^aOS, where Ccai is 0.05 and 0.12 for the MIPS 70 fim ([Gordon et al][2007| 



and MIPS 160 /im ( Stansberry et aL]|2007 ) bands respectively. The few sources (9) requiring large 
aperture measurements have additional measurement errors of 10-20%. 

We find that the positional fitting errors calculated by APEX are 2.0 times larger on average 
than the expected radial positional errors of ~ 0.6(FWHM/SNR) given by Condon (19971 for 



all SNRs. We adopt the APEX errors, but treat them as 2a uncertainties. For these data the 



relationship of Condon (19971 is appropriate for typical SNRs, but underestimates the errors for 
the highest SNRs. We adopt a floor on the uncertainty of epos = 0.5 pixel (~ 0.1 FWHM, Table 2), 



representing the difficulty in deriving positions to better than a fractional pixel regardless of the 
SNR. The cataloged radial positional errors {2a uncertainties) are given by (o"^ -|- a^ 
where ax and ay are the fitting errors in the x-y plane derived by APEX. 
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5. Results and Discussion 
5.1. Source Counts 

The source counts were derived from SNR> 5.0 source lists corrected for completeness over the 
region within the images with high coverage (> 100 at 70 //m and > 20 at 160 /^m), corresponding 
to an effective survey area for the 70 and 160 fields of 2.471 and 2.438 sq-deg, respectively. 
The public catalogs presented in Sec. 4.2 include sources outside of this nominal area with lower 
coverage. The effective completeness value for each flux density bin was computed by integrating 
over the completeness curves (Fig. 3&4) as a function of flux density across each bin with weights 
based on the measured slope of the source counts. 

Figures 5&6 show the Euclidean-normalized differential source counts {dN/dS x 5^'^) at 70 /xm 
and 160 /um, and the results are tabulated in Tables 6&7, respectively. The counts are calculated 
for independent flux density bins. The error bars represent the Poisson errors associated with the 
number of sources and the uncertainties on the completeness corrections. The grey-region shows 
the range of values implied from the most statistically accurate counts within each flux density 
bin, including the additional uncertainty due to the systematics associated with the calibration of 
MIPS (5% at 70 /xm and 12% at 160 /xm). 

At 70 /um the S-COSMOS counts are measured down to a level of 10 mJy which is near the peak 



of the Euclidean-normalized differential source counts (Frayer et al. 2006b). The faint (~ lOmJy) 



S-COSMOS counts at 70 //m (Fig. 5) agree with those measured for GOODS-N (Frayer et al 



2006bD and the xFLS ( [Frayer et aL]|2006a[ ). The bright counts for S-COSMOS also agree with 
the measurements from the SWIRE survey. At intermediate flux densities of around 20-30 mJy, 
the counts for S-COSMOS are slightly lower than those found for the xFLS and model predictions 



of Lagache et al. (2004 1 . All the counts have been placed on the same scale by matching the 
calibration, color corrections, and the PRE adopted for S-COSMOS. The correction factors for the 
other data sets are given in Table 5. The previous results for the xELS, SWIRE, and GOODS-N 
did not include the aperture correction for the flux density outside of the measured PRF. The 
counts for SWIRE are based on the public catalogs (2005 November, Data Release 3 [DR3] ) which 
cover 49 sq-deg. The SWIRE counts at 70 or 160 /xm have not previously been published and are 
presented here to constrain the counts at the brightest flux densities. Only the high SNR (> 10) 
SWIRE sources with completeness levels near one are presented in Figures 5&6. 

At 160 /xm the S-COSMOS counts are measured down to a level of 60mJy. Measurements of 
deeper counts are limited by confusion (Sec. 5.2). The measured counts at 160 /xm agree with the 
counts measured previously in the xFLS ( [Frayer et al. 2006a) and the counts derived here based 



on the SWIRE survey. As done for 70 /um, all counts are placed on the current calibration scale 
(Table 5). These correction factors also include a decrease in the flux densities due to the new 
PRF derived here (Sec. 3.4). The previous PRF images used by the xFLS and SWIRE surveys are 
affected by the flux non-linearity for bright sources at 160 fim. 
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The faintest (5*160 < 80mJy) and brightest (S'160 > 500 mJy) source counts are consistent 

model, but the observed counts for all three surveys at intermediate 
flux densities are about a factor of 1.5 times lower than the model implies. The observed 160 
/im counts show a steep increase in the differential counts for decreasing flux densities {dN/dS oc 
5.-3.5±o.2) for 5160 < 150 mJy. The slope for the faint (10-20 mJy) S-COSMOS 70 /im source 
counts {dN/dS oc S*"^'^^"'^) is roughly consistent with the 160 /xm slope and is slightly steeper 
than the model predictions. The differences between the model and the source counts at 70 ^va. 
is not as clear as those seen at 160 /xm. The results for 160 /xm may suggest the importance of 
galaxies having more cold dust than assumed in the models. Observations with the future Herschel 
telescope should provide better constraints on the FIR SEDs of distant galaxies. 



with the Lagache et al. (|2004 



5.2. Confusion Level of the MIPS 70 and 160 /im Bands 



The S-COSMOS 70 /im data are dominated by instrument uncertainties and are not deep 



enough to constrain the confusion noise. Frayer et al. (2006b) measured the MIPS 70 /tm confusion 



noise based on the much deeper 70 /im observations of GOODS-North. Based on [Frayer et al. 

tors (Table 5), th( 
0.35 ±0.15 mJy {q 



(2006b I and the calibration scaling factors (Table 5), the updated extragalactic confusion noise 

5), including the updated systematic 



level for the MIPS 70 /im band is a, 
error on the flux calibration. 

In contrast to the 70 /im data, the S-COSMOS 160 /tm data contain a significant noise compo- 
nent due to confusion. We estimate the confusion noise at 160 /xm following the empirical technique 



performed at 70 /im (Frayer et al. 2006b). A direct empirical measurement of the confusion noise 



for the MIPS 160 /im band has not been published previously. Dole et al. (2004) report a confusion 



level at 160 /im based on the techniques of Dole et al. (2003) and the models of Lagache et al. (2004) 



They find that the predicted confusion noise from their models is in reasonable agreement with the 
observations, assuming that the instrumental noise (ci/) follows the theoretical MIPS model and 
integrates down as ct/ oc t~^'^. However, the early processing of MIPS-Ge data was not optimal, 
and an accurate empirical measurement of the total instrumental noise (including photon noise, 
detector noise, and noise associated with the data processing) is required to measure the confusion 
level. 

The instrument noise was estimated empirically by subtracting pairs of data subsets with the 
same integration time and covering the exact same region on the sky (which removes sources and 
any remaining cirrus structure after filtering). We fit the instrument noise measurements as a 
function of integration time for combinations of deep pairs of data sets and obtain aj oc i-0-49±o.03_ 
This result is consistent with idealized data {a oc t~^'^) and highlights the success of the reduction 
methods in removing systematic artifacts. The extrapolation of instrument noise for half of the 
data to the full data set yields ci/ = 0.1134ib0.0033MJysr~^, where the uncertainty represents the 
rms measurement error combined in quadrature with the error associated with the extrapolation. 
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Following the terminology of Dole et al. (20031, the total noise {ct) represents the noise after 
the extraction of sources above a limiting flux density {Sum), and the photometric confusion noise 
(cTc) represents fluctuations due to sources with flux densities below Sum- The confusion noise is 
given by ac = {cr^ — aj)^'^, which is appropriate for an approximate Gaussian distribution of the 
noise after the extraction of sources. We iterate between source extraction at different limiting flux 
densities and confusion noise measurements until we converge to a solution with q = Sn^/ac = 5. 
For the g = 5 solution, we derive ax = 0.1772±0.0089MJysr-i and ac = 0.1362±0.0119MJysr"\ 
for a limiting source flux density of 5160 = 50mJy. 

Noise measurements in surface brightness units (MJyrsr~^) depend on the pixel scale, and 
all measurements here are based on the 8" pixel scale of the 160 fim image. The uncertainties 
on the noise measurements given in MJysr"^ also do not include the 12% systematic calibration 
uncertainty. Including the systematic calibration uncertainty and the measured conversion between 
surface-brightness noise and point-source noise (Sec. 4.1), the point-source confusion level of the 
MIPS 160 fxm band within the COSMOS field is ac = 10.0 it 3.1 mJy. The systematic uncertainties 
contribute 71% to the total error budget, while the random errors contribute 29% to the total error 
budget. 

The measured confusion noise is not entirely due to galaxies. Unlike the case for the MIPS 70 



Dole et al. 


2003 


Frayer et al. 


2006b 



We estimate the level of confusion due to Galactic cirrus for the MIPS 160 iim band using the 
background estimates from the Spitzer tool for planning observations (SPOT) which is based on 



the interstellar medium (ISM) maps from the Diffuse Infrared Background Experiment (Schlegel et 
al] |1998[ ). At the effective wavelength of 155.9 fim of the MIPS 160 j-im band, the ISM background 
in the direction of the COSMOS field is about 2MJysr~^. Using the calculations of Jeong et al. 
(2005), this background level corresponds to an ISM confusion noise of (Tc(ISM) ~ 3.4 mJy. Hence, 
the confusion noise due to unresolved galaxies is (Tc(gal) = ((T^(total) — cr^(ISM))°-^ = 9.4ib3.3mJy. 



The derived confusion limit agrees fairly well with predictions. Dole et al.| (2004 1 report a 
Source Density Criterion (SDC) limit of 40mJy {adqsDC = 3.8) = 10.6 mJy) based on the model 



predictions of Lagache et al. (2004). For a direct comparison with the empirical measurement, the 
Dole & Lagache et al. predictions suggest adq = 5) = 12.5 mJy. We find an extragalactic confusion 
level for the MIPS 160 ^m band of ac{q = 5) = 9.4 it 3.3 mJy, which is just slightly lower than the 
predicted value. 



6. Concluding Remarks 

We present the 70 and 160 iJ,m observations of the COSMOS field and describe the products. 
This is the first extragalactic survey available to the public at 70 and 160 /im that has been placed 
on the calibration scale derived from the recent MIPS calibration papers. We provide updated 
correction factors for the previously released catalogs of the xFLS ( Frayer et al.|2006a ) and SWIRE 
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(Lonsdale et al. 2004) programs. Counts are presented based on the S-COSMOS and previous 



surveys and are found to be in reasonable agreement with the model of Lagache et al. (2004 1 . 
However, the faint 160 //m source counts are significantly steeper than model predictions. We 
measure an empirical extragalactic confusion noise level of cTc = 9.4 it 3.3 mJy {q = 5) for the 
MIPS 160 /im band. In comparison, the expected confusion noise at 160 ^iui for Herschel is about 
1.5-2 mJy {Herschel Observation Planning Tool, HSPOT, version 3.4). Future observations with 
the Herschel telescope should constrain the counts and far-infrared properties better than can be 
done currently due to confusion for the MIPS 160 /im band. 

We thank our colleagues associated with the Spitzer mission who have made these observations 
possible. This work is based on observations made with the Spitzer Space Telescope, which is 
operated by the Jet Propulsion Laboratory, California Institute of Technology under a contract with 
NASA. Support for this work was provided by NASA through an award issued by JPL/Caltech. 
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Table 1. Spitzer COSMOS Observations with MIPS 



Epoch 


Dates 


Campaign 


AORs 


Cycle-2 


2006 Jan 06-09 


MIPS006300 


21 


Cycle-3a 


2007 Jan 01-13 


MIPS010800 


70 


Cycle-3b 


2007 May 18-28 


MIPS011900 


61 


Cycle-3c 


2008 Jan 06-07 


MIPSO 13500 


5 



Note. — Spitzer programs 20070 and 30143. 



Table 2. Summary of Product Properties 



Property 70 /nm 160 /xm 



Image pixel scale 4" 8" 

Image FWHM resolution 18.6" 39" 

Nominal coverage level 100 20 

Median coverage" 140.3 28.3 

Median exposure time" 1350 s 273 s 

Surface-brightness noise (la)" 0.13 MJysr"^ 0.18 MJysr"^ 

Point-source noise (Icr)" 1.7 mJy 13mJy 

Area mapped" 2.471 sq-deg 2.438 sq-deg 

Number of sources (> S.Oa) 1512 499 



Note. — "Measured properties for coverages (number of obser- 
vations) larger than the nominal coverage level. 
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Table 3. Spitzer COSMOS 70 /xm Catalog 



Source Name 


a(J2000) 


J(J2000) 


Error 


570 


a{S70) 


SNR 


Flag 




(deg) 


(deg) 


(arcsec) 


(mJy) 


(mJy) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


SCOSMOS70 J095546.9+013605 


148.945497 


1.601560 


5.0 


20.6 


4.9 


5.3 


P 


SCOSMOS70 J095552.5+014243 


148.968967 


1.712061 


2.2 


226.2 


21.9 


21.3 


P 


SCOSMOS70 J095605.9+014443 


149.024976 


1.745507 


4.7 


17.1 


3.9 


5.7 


P 


SCOSMOS70 J095610.1+014213 


149.042376 


1.703652 


4.8 


14.3 


3.4 


5.4 


P 


SCOSMOS70 J095615.0+014315 


149.062575 


1.720973 


4.0 


16.6 


3.2 


7.1 


P 



Note. — The entire S-COSMOS 70 fim catalog (> 5. Oct) is presented in the electronic edition of the 
Astronomical Journal. A portion of Table 3 is shown here for guidance regarding the form and content of 
the catalog. Column (1) gives the source names following the lAU designations. Column (2) and column 

(3) are the right ascension (a) and declination (S) J2000.0 source positions in decimal degrees. Column 

(4) is the radial positional uncertainty in arcsec (2ct). Column (5) is the total flux density measurement 
in mJy. Column (6) is the flux density error in mJy (lo"), including the systematic uncertainty of the 
absolute flux density scale. Column (7) is the SNR of the peak. Column (8) is the flag for the flux density 
measurement method. A flag of "p" indicates that the source was fitted by the PRF, "a" indicates an 
aperture measurement with a diameter of 96", and "ap" indicates an aperture measurement whose flux was 
divided between two components based on the relative strengths of their peaks. 



Table 4. Spitzer COSMOS 160 /im Catalog 



Source Name 


a(J2000) 


(5(J2000) 


Error 


5160 


(t(5160) 


SNR 


Flag 




(deg) 


(deg) 


(arcsec) 


(mJy) 


(mJy) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


SCOSMOS160 J095619.9+012425 


149.083255 


1.407074 


9.2 


134.9 


40.1 


5.6 


P 


SCOSMOS160 J095621.2+015951 


149.088411 


1.997640 


9.7 


122.5 


37.2 


5.4 


P 


SCOSMOS160 J095635.5+014332 


149.147947 


1.725679 


9.2 


122.7 


35.4 


5.9 


P 


SCOSMOS160 J095635.7+012544 


149.149042 


1.429048 


4.7 


1073.7 


188.5 


18.0 


P 


SCOSMOS160 J095636.9+015415 


149.153962 


1.903362 


5.0 


1093.4 


209.8 


13.9 


P 



Note. — The entire S-COSMOS 160 /xm catalog (> 5. Oct) is presented in the electronic edition of the 
Astronomical Journal. A portion of Table 4 is shown here for guidance regarding the form and content of the 
catalog. Column headers arc the same as described in Table 3 for columns (l)-(7). Column (8) is the flag for 
the flux density measurement method where "p" indicates a PRF measurement and "a" indicates an aperture 
measurement with a diameter of 4'. 
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Table 5. Scaling Factors for MIPS-Ge Surveys 



Survey, Version 



70 fim 160 i-im 



S-COSMOS, version 1&2 
S-COSMOS, version 3 



xFLS(Frayer et al. 



2006a I 



SWIRE (2 005 November, DR3 ) 
GOODS-N (Frayer et al.|[2006b| ) 



1.25 
1.00 
1.20 
1.10 
1.15 



1.07 
1.00 
0.97 
0.98 



Note. — Multiplicative scaling factors required to 
place catalogs on the same calibration scale derived 



for the official calibration of MIPS ( [Gordon et al. 



2007; Stansberry et al. 20071 and used here for the 



S-COSMOS, version 3 products. Factors include the 
combination of color-corrections, calibration updates, 
PRF updates, and the aperture correction for emis- 
sion outside of the measured PRF. 
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Table 6. S-COSMOS TO/xm Source Counts 



Observed 






Observed 




Eddington 




Average 


Slow 


Shigh 


Number 


Completeness 


Bias 


dN/dSS"^-^ 


(mJy) 


(mJy) 


(mJy) 








(gal sr~^ Jy^'^) 


11.02 


10 


12 


245 


0.68± 0.05 


1.06 


2872± 280 


12.84 


12 


14 


186 


0.90± 0.03 


1.04 


2478± 199 


15.30 


14 


17 


180 


0.98± 0.02 


1.03 


2304± 178 


19.08 


17 


22 


144 


l.OOib 0.02 


1.01 


1937± 166 


25.76 


22 


30 


108 


l.OOib 0.02 


1.01 


1924± 189 


36.26 


30 


45 


92 


1.00± 0.02 


1.00 


2086± 221 


54.05 


45 


70 


46 


l.OOi 0.02 


1.00 


1697± 252 


86.36 


70 


110 


22 


l.OOi 0.02 


1.00 


1637± 350 



Note. — The counts are based on the > 5. Oct catalog and were measured for an area of 
2.471 sq-deg with a coverage of greater than 100 observations. The counts were corrected 
for completeness and the Eddington bias (5'obs/'S'true)- The tabulated uncertainties for the 
differential counts {dN/dS S'^-^) include the Poisson noise and the completeness uncertainties, 
but do not include the additional systematic uncertainty (5%) associated with the calibration 
of the MIPS 70 /xm band. 
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Table 7. S-COSMOS 160 fim. Source Counts 



Observed 






Observed 




Eddington 




Average Si, 


Slow 


Shigh 


Number 


Completeness 


Bias 


dN/dSS^-^ 


(mJy) 


(mJy) 


(mJy) 








(gal sr~^ Jy^'^) 


65.05 


60 


70 


82 


0.63± 0.05 


1.09 


16619±2260 


74.82 


70 


80 


83 


0.85± 0.03 


1.08 


17875±2060 


89.01 


80 


100 


83 


0.95± 0.02 


1.06 


12713±1420 


139.18 


120 


160 


40 


l.OOi 0.02 


1.02 


9445±1505 


187.72 


160 


250 


29 


l.OOi 0.02 


1.01 


6526±1218 



Note. — The counts are based on the > 5. Oct catalog and were measured for an area 
of 2.438 sq-deg with a coverage of greater than 20 observations. The counts were corrected 
for completeness and the Eddington bias (S'obs/'S'true)- The uncertainties for the differential 
counts {dN / dS S'^-^) include the Poisson noise and the completeness uncertainties, but do 
not include the additional systematic uncertainty (12%) associated with the calibration of 
the MIPS 160 fiia band. 
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151.0 150,5 150.0 149.5 149.0 
RA (J 2000.0) 

Fig. 1. — The S-COSMOS 70 /7,m image shown on a logarithmic grey-scale. The solid lines show 
contours representing coverages of 50, 100, and 250 observations. The central rectangle (highest 
coverage > 250) shows the Cycle-2 test field, and the dotted line shows the extent of the RST ACS 
field. 
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151.5 151.0 150.5 150.0 149.5 149.0 

RA (42000.0) 

Fig. 2. — The S-COSMOS 160 /xm image shown on a logarithmic grey-scale. The solid lines show 
contours representing a coverage of 10, 20, and 50 observations. The deep central contour shows 
the Cycle-2 test field, and the dotted line shows the extent of the HST ACS field. 
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Fig. 3. — The 70 /7,m completeness levels (SNR > 5.0) as a function of flux density and coverage. 
The effective average completeness level used for the source counts is given by the solid line (Total 
Cov> ICQ). Simulation results for the deep area (Gov > 260) and typical coverage values (130 < 
Gov < 150) are shown by the grey diamonds and squares. The dashed and dotted lines represent 
the expected completeness level based on the SNR threshold for the regions used in the simulations. 
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50 100 150 200 

Point Source Flux Density [mJy] 



Fig. 4. — The 160 /Ltm completeness levels (SNR > 5.0) as a function of flux density and coverage. 
The effective average completeness level used for the source counts is given by the solid line (Total 
Gov > 20). Simulation results for the deep area (Gov > 50) and typical coverage values (25 < 
Gov < 33) are shown by the grey diamonds and squares, along with the expected completeness 
level based on the SNR threshold and coverage values (dashed and dotted lines respectively) 
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Fig. 5. — The Euclidean-normalized differential source counts at 70 fim compared to the evolu- 
tionary model of Lagache et al. (2004 1 (solid line). The S-COSMOS and SWIRE counts are new 
measurements. The xFLS and GOODS-N counts are from Frayer et al. ( 2006a|b | respectively. All 
data are on the same calibration scale. The grey region shows the range of values for the accurate 
data points within each flux bin and includes the additional systematic uncertainty associated with 
the calibration at 70 /^m (5%). 
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Fig. 6. — The Euclidean-normalized differential source counts at 160 fim compared to the evolu- 
tionary model of Lagache et al. (2004 1 (solid line). The S-COSMOS and SWIRE counts are new 
measurements, while the xFLS counts are from Frayer et al. (2006a I . All data are on the same 



calibration scale. The grey region shows the range of values for the accurate data points within 
each flux bin and includes the additional systematic uncertainty associated with the calibration at 
160 (12%). 



